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The superconducting properties of carbon nanotube ropes are studied using a new computational 
framework that incorporates the renormalization of intratube interactions and the effect of intertube 
Coulomb screening. This method allows to study both the limits of thin and thick ropes ranging from 
purely one-dimensional physics to the setting of three-dimensional Cooper-pair coherence, providing 
good estimates of the critical temperature as a function of the rope physical parameters. We discuss 
the connection of our results with recent experiments. 
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Introduction. — Carbon nanotubes offer nowadays a 
great potential for technological applications. They dis- 
play semiconducting or metallic properties depending on 
the helicity of the tubule [1] , what is a useful feature to 
build electronic devices. The electron correlations are 
significant in the metallic nanotubes [2-5], explaining 
the experimental observation of features characteristic of 
Luttinger liquid behavior [6,7]. Quite remarkably, super- 
conducting correlations have been also observed, in the 
form of the proximity effect [8,9] as well as in the sharp 
decrease of the resistance at low temperatures [10]. 

The superconducting transitions reported in Rcf. [10] 
have been measured in various samples of nanotube 
ropes. These are supposed to be made of a large num- 
ber (~ 300) of single-walled nanotubes, and in all cases 
the critical temperature has been found below 1 K. In a 
different kind of experiment [11], strong superconducting 
correlations have been observed in individual nanotubes 
of very small diameter placed in a zeolite matrix. It has 
been claimed that such correlations should correspond to 
a critical temperature of the order of ^ 15 K [11]. 

An explanation of the superconductivity (SC) in the 
nanotube ropes has been already proposed [12,13]. It has 
been argued that, as long as the Cooper pairs are formed 
at zero total momentum, they do not find the obstacle 
that single electrons have to tunnel between neighbor- 
ing nanotubes due to the misalignement of the lattices 
[14,15]. The existence of a superconducting transition is 
possible in the ropes as the coherence in the transverse 
directions is established through the tunneling of Cooper 
pairs between the nanotubes [13]. In this scenario it has 
been assumed that the screening of the Coulomb repul- 
sion is the dominant effect in large nanotube ropes and 
that the backscattering and Umklapp interactions are not 
essential to explain the SC of these systems. 

In this Letter we propose a new computational frame- 
work to deal simultaneously with the Coulomb screening 
due to the intertube interactions and the scaling of all the 
intratube interactions in the low-energy theory. In this 
way we will obtain a precise estimate of the critical tem- 
peratures that can be reached. The interplay between 



all the possible interactions in the scaling procedure will 
allow us to investigate the development of 3D coherence 
in the rope, as well as to discover new phases which cor- 
respond to different charge instabilities in the system. 
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FIG. 1. Plot of the low-energy linear branches of carbon 
nanotubes. The spinors in the inside (outside) correspond 
to the relative electron amplitudes in the two sublattices of 
graphene rolled up to form a zigzag (armchair) nanotube. 

Electron-phonon interactions. — The effective in- 
teraction coming from the exchange of phonons may lead 
to attraction or repulsion between the electrons, depend- 
ing on the particular channel of interaction considered 
[16,17]. We will follow the convention of classifying the 
four-fermion interactions into different channels with re- 
spective coupling constants g'f'^ [18]. The lower index 
discerns whether the interacting particles shift from one 
Fermi point to the other {i — 1), remain at different 
Fermi points (1 = 2), or they interact near the same Fermi 
point {i = 4). The upper label follows the same rule 
to classify the different combinations of left-movers and 
right-movers, including the possibility of having Umk- 
lapp processes (j = 3). 

Let us focus on the opposite cases of armchair and 
metallic zigzag nanotubes [1]. All of them have two pairs 
of linear branches in their band structure, crossing at 
respective points at the Fermi level (for undoped nan- 
otubes) as shown in Fig. 1. Obviating the momentum de- 
pendence of the electron modes, the relative amplitudes 
between points in the two sublattices of the honeycomb 
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lattice arc given by the spinors in Fig. 1 [19]. 

Using the fact that the electron modes with opposite 
momenta have complex conjugate amplitudes, it can be 
shown that the contributions from phonon exchange to 
9'\\ 92^^ ^iid ^4^^ are negative, for frequencies 
below the characteristic phonon energy. However, the 
contributions from Umklapp processes to g^\g'^'' and 

(3) 

94 have the opposite sign, since the respective products 
of electron-phonon couplings turn out to be negative [17]. 
We recall that the electron-phonon couplings are given by 
a sum involving products of electron modes in ncarest- 
neighbor sites [16]. Thus, the contributions to g\ and 

g'^'p have attractive character in the zigzag nanotubes, 
due to the relation of conjugation between the electron 
modes, but they have repulsive character in the armchair 
nanotubes, since only one of the two scattered modes 
reverses its sign there when shifting to a neighboring site. 

The signs of the different contributions are based on 
symmetry rules, and they hold irrespective of whether the 
exchanged phonons are optical or acoustic. The energy 
of the latter at momentum 2k p is close to the energy of 
optical phonons near fc = 0, so that we will assume equal 
strength in absolute value, \g\, for all the four-fermion 
interactions arising from the exchange of phonons with a 
typical Debye frequency between 0.1 and 0.2 eV. 

We will deal with effective interactions coming from 
the exchange of phonons within the same nanotube. In- 
teractions mediated by the exchange of phonons in differ- 
ent nanotubes require mainly the coupling of out-of-plane 
phonons, which have a much lower energy scale. These 
effects have been studied in Ref. [20], with the result that 
the intratube superconducting order parameter is more 
relevant than the intcrtube one, when the system is in 
the regime with dominant attractive interactions. 

Coulomb screening. Regarding the Coulomb in- 
teraction, it has been shown that the strength of the 
backscattering and Umklapp processes mediated by the 
Coulomb potential is reduced by a relative factor 
0.1 a/i?, in terms of the ratio of the lattice spacing a to 
the nanotube radius R [2-4]. This means that, for typ- 
ical nanotubes in a rope, these processes can be safely 
neglected in favor of the contributions from phonon ex- 
change, whenever the strength in absolute value [^l of 
the latter is above ~ O.OSuf. Therefore, the competition 
between the long-range Coulomb interaction and the ef- 
fective interaction from phonon exchange takes place in 
the channels with small momentum transfer correspond- 
ing to the couplings \ S'2^\ fitf ^ and gf^. 

Dealing for the time being with such small momentum- 
transfer interactions, the dressed vertices are given by 
self-consistent equations of the type shown in Fig. 2. 
It has to be realized that, in the case of a rope, there 
are interaction vertices o'f ^ for currents within the same 
nanotube as well as interaction vertices K,^"'' between cur- 



rents at different metallic nanotubes. If the number of 
them is n, this is also the number of different terms in 
the sum that appears in Fig. 2. 





FIG. 2. Diagrams contributing to the screening of Coulomb 
interactions between currents la, h of well-defined chirality 
and spin inside a rope. The second term takes into account 
the polarization of the n metallic nanotubes in the rope. 

The self-consistent equations for the intratube 

(i) 

and the intertube V-^ ' vertices can be solved in the 

same fashion as for the simpler case of the Luttinger 
model [21]. In order to obtain a solution in the form 
of coupling constants, we have taken the static approx- 
imation in the polarization operators. We have also 
chosen an average strength for the Coulomb potential 

V — (e^/27r) log |/cc/fcol) where kc is the short-distance 
cutoff and ko ~ 10~^A;c accounts for the long-distance 
average. 

The appropriate bare values of the couplings are g^^ = 

(2) (2) , 1 (4) 1 

V + g , gl = V g , = + .9 .92 = v =p g, where 
g < is the contribution from phonon exchange. The 
upper sign applies to the armchair nanotubes and the 
lower sign to the zigzag nanotubes. In the former case, 
the solution of the self-consistent equations is 
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In the case of the zigzag nanotubes, we have 
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for all of the four dressed couplings. The effect of the 
Coulomb screening is well captured in the expression of 
the dressed vertices. The RPA-like diagrams considered 
screen repulsive interactions and enhance attractive in- 
teractions, and the latter always prevail for sufficiently 
large values of n. For Algl/eoT^vp close to 1, strong non- 
perturbative effects arise from the same kind of singular- 
ity already observed in the bosonization picture [22,13]. 
The limit n — > 0, 5' — » recovers the bare vertices used 
in precedent studies of SC in isolated nanotubes [17,18]. 

Intratube scaling equations. — The preceding 
physical description has to be completed by incorporat- 
ing the effect of the backscattering and Umklapp effective 
interactions arising from phonon exchange. We set the 
signs of their respective initial strengths (equal to \g\ in 
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absolute value) according to the above discussion, and 
take into account their nontrivial scaling behavior when 
the frequency probe is progressively reduced. The cou- 
plings scale according to anomalous dimensions which 
can be written in terms of the dressed couplings D^-'K 
To second order in the backscattering and Umklapp cou- 
plings, the scaling equations read 
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where / stands for (minus) the logarithm of the energy 
scale measured in units of the high-energy cutoff Ec (of 
the order of ~ 0.1 eV). These equations are similar to 

those obtained in Rcf. [18], except for the replacement 



of ^2^^ , 5(2^ ■* , g'f'' and g)^' by the respective dressed cou- 
plings. We apply this replacement to the rest of the 
equations in the model, including those for the small 
momentum-transfer couplings, which we do not write 
since they are similar to those in Ref. [18]. In this way we 
incorporate the finite renormalizations arising from the 
interaction among the metallic nanotubes in the rope, 
which set the initial values of the dressed couplings in 
accordance to Eqs. (l)-(3). 

Intertube coherence. — The ropes are made of nan- 
otubes with different helicities and diameters, and this 
compositional disorder frustrates the formation of any 
phase with spin-density-wave or charge-density-wave or- 
der in the rope. On the other hand, the nanotubes 
may have strong superconducting correlations. Their 
strength is given by the combination of the couplings 
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and g\ , which govern the propagation 
of the Cooper pairs along each nanotube. The interac- 
tions corresponding to 5^^^ and g^^ are always attractive, 
pointing at the enhancement of the superconducting cor- 
relations with s-wave symmetry at low energies. 

In the compositionally disordered ropes, the coherence 
is established by the tunneling of Cooper pairs between 
neighboring metallic nanotubes. The amplitude J for 
this process has a dependence on the energy scale char- 
acterized by the anomalous dimension A, given by the 
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combination D.^ 
tion for J is [23] 



dJ/dl = AJ + cAI 



The scaling equa- 
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where It is the intertube hopping for each particle of 
the Cooper pair, of the order of ~ 0.01 eV. We re- 
call that the tunneling amplitude for single electrons is 
in general much smaller than this quantity, due to the 
mismatch between the Fermi points of neighboring nan- 
otubes [14]. This introduces a factor of suppression that 
sets the single-particle intertube hopping about three or- 
ders of magnitude below the estimate for tr [15]. 

If the generalized backscattering and Umklapp pro- 
cesses are fixed to zero, the system is at the fixed point de- 
scribed in the bosonization picture [13]. That fixed point 
is in general unstable against chirality-brcaking pertur- 
bations by which the system flows to a strong-coupling 
fixed point. Nevertheless, the crossover behavior is very 
slow, being characterized by a relatively small crossover 
frequency w^o- In a large region of the phase diagram the 
flow is cut off by the onset of the 3D coherence before w^o 
has been reached. In these conditions, the SC phase re- 
mains in the region of influence of the weak-coupling fixed 
point, which governs then the main physical properties. 

Results and discussion. Following the above ap- 
proach, we characterize the superconducting transition 
by looking for a pole in the 3D propagator 'R,{uj) of 
the Cooper pairs along the rope. This is related to the SC 
response function T^*^"' (w) in the individual nanotubes by 
the Schwinger- Dyson equation 1/TZ{lj) = 1/TZ^^\lj) — J, 
where J is the (dimensionless) pair hopping measured in 
units of the high-energy cutoff. 

Moreover, we can study the competition between the 
SC and other charge instabilities in the rope by re- 
bosonizing the system at each step in the integration of 
Eqs. (4)- (8), considering backscattering and Umklapp 
couplings as small perturbations. The generalized Lut- 
tinger liquid parameters of the individual nanotubes take 
the form 



nvF + {Dt^ ± DD - iD2' ± D4') 
\nvF + {Di'^±Di'y) + {D^^^±D^:^) 
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We obtain complementary information about the physi- 
cal properties of the system by studying the regions where 
the K± parameters either vanish or diverge. This hap- 
pens when some of the functions given by Eq. (10) de- 
velops a branch-cut at some energy scale, what has to be 
understood as the onset of a phase transition in the sys- 
tem. The different phases which arise in this approach 
have been represented in Fig. 3, where contour lines of 
constant critical scale Wc have been also plotted. 

The weak-coupling character of the SC phase allows to 
estimate the transition temperatures from the values of 
the pole cOc shown in Fig. 3. The comparison with the ex- 
periments reported in Ref. [10] can be made by recalling 
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that the ropes used there have about 100 metalhc nan- 
otubes, while the transition temperatures measured are 
always below 1 K. When compared to the high-energy 
cutoff, that value corresponds to Ic ~ 7. From the re- 
sults shown in Fig. 3, we find that the strength of the 
attractive interaction 4:\g\/nvF has to be « 0.3. 
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FIG. 3. Phase diagram of undoped armchair (left) and 
doped zigzag (right) nanotubes as a function of the number 
of metallic nanotubes n in the rope and the strength of the 
attractive interaction G — Algl/TTVp- Thick lines represent 
the phase boundaries between a 3D phase-coherent supercon- 
ducting phase (SC) and the ID phases characterized by the 
breaking of the parameters and K- (the narrow phase to 
the left of the second diagram corresponds to the divergence 
of the latter). Thin lines are contours of constant critical fre- 
quency LOc, starting from above with \n{Ec/uJc) = 2, 4, 6, . . . 

Given the above estimate of g, we conclude that the 
superconducting correlations have to disappear below a 
number of metallic nanotubes n < 50, in ropes of the 
type measured in Ref. [10]. That is, for nanotubes with a 
typical diameter of w 1.4 nm, the effect of the Coulomb 
interaction overcomes any source of effective attraction 
for that range of n. This is in agreement with the experi- 
mental measures carried out by the collaboration of Ref. 
[10] in thin ropes with about 40 nanotubes [24]. 

The lower boundary of the SC phase in Fig. 3 almost 
coincides with the boundary at which the pair-hopping 
amplitude ceases to grow large at low energies. When this 
happens, the nanotubes in the rope start to behave as un- 
coupled ID systems, what explains the disappearance of 
the superconducting transition. A phase with vanishing 
compressibility appears for low values of [^l in the un- 
doped system. The fact that the regime with K+ = is 
only present at half-filling makes clear that it corresponds 
to the extension of the Mott insulating phase found in the 
studies of carbon nanotubes with respulsive interactions 
[4,5]. In the case of doped nanotubes, the repulsive Umk- 
lapp interactions are not present in the system, and the 
decoupling of the nanotubes in the rope leads instead to 
a phase with divergent compressibility in the channel of 
the total charge, as shown in Fig. 3. 

Another phase with divergent compressibility arises at 



strong coupling in the channel of the mismatch of charge 
in the two gapless subbands. These divergences have 
the same character of the Wentzel-Bardeen singularity, 
which is the natural counterpart of the superconducting 
instability [22,20]. A divergent compressibility implies 
that the density-density correlations become increasingly 
large, pointing at the onset of a regime where there is no 
homogeneous descripition of the system. For sufRciently 
strong attraction, the electronic charge finds then more 
favorable its macroscopic segregation than the formation 
of Cooper pairs. We have shown that such a phase may 
be observable in ropes, either as a consequence of a strong 
attractive interaction or, in doped systems, from the de- 
coupled behavior of the metallic nanotubes at weak cou- 
pling, for temperatures slightly below 1 K. 
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